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Abstract 

The nearest-neighbor mass-spacing distribution of the meson and baryon spectrum (up to 2.5 
GeV) is described by the Wigner surmise corresponding to the statistics of the Gaussian orthogonal 
ensemble of random matrix theory. This can be viewed as a manifestation of quantum chaos in 
hadrons. 
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Quarks are glued into colorless objects called hadrons, forming a large and complex 
hadron spectrum. This spectrum is the main source of information about the quark-gluon 
interactions in the confinement regime, and hadron spectroscopy has been receiving a lot of 
attention both theoretically and experimentally. 

In the study of complex spectra in general, statistical methods proved to be useful. One 
of the early successes of statistical analysis was Wigner's celebrated discovery of the fact 
that fluctuation properties of complex nuclear spectra are described by the Random Matrix 
Theory (RMT) ||, |, ^, ||. It was later realized that these properties, referred to as 
Wigner-Dyson properties, have a great deal of universality and appear in spectra of many 
physical systems, ranging from quantum dots to lattice gauge theories [0, H, ^ p!0| . 

In this Letter we show that the experimentally measured hadron spectrum has at least 
some Wigner-Dyson properties as well. Once this is established the link to the underlying 



dynamics is made via the Bohigas-Giannoni-Schmit (BGS) conjecture [TTp, which states 
that the Wigner-Dyson properties are generic to spectra of the systems with a quantum 
analog of chaotic dynamics — the "quantum chaos" (see, e.g., [|1^, 0). 

Recent theoretical work in lattice QCD 0, § ^ and quantum-mechanical Yang-Mills 
models [|1^, 0] has already addressed the statistical properties of QCD spectrum as well 
as the quantum chaos aspects. Our findings will provide certain empirical support of these 
results and, hopefully, motivate further studies in this direction. 

We look at the experimentally measured mass spectrum of hadrons up to 2.5 GeV taken 
from the Particle Data Group (PDG) Summary Tables More specifically, we consider 
A^, A, A, and S baryons up to A^(2200), and all the mesons listed in the Summary Tables 
up to /2(2340). In doing so we exclude the poorly known "one-star baryons." We have only 
verified that they do not have a significant impact on our results. 

The spectrum can be organized into multiplets, see Fig. |l], characterized by a set of definite 
quantum numbers (QNs): isospin, spin, parity, strangeness, baryon number, and, in the case 
of mesons, charge conjugation. Using these data we would like to examine the probability 
distribution of spacing (mass-splitting). Si = rrii^i — nii, between nearest-neighbor hadrons 
within one multiplet (same QNs). We thus do not need to consider channels with only a 
single state below 2.5 GeV. 

Each multiplet provides a number of spacings which are gathered into a common array 
of spacings. We actually will be considering three arrays: {Si}B, {Si}M, and {Si}H = 
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FIG. 1: Light hadron spectrum (excluding 1-star baryons) in the form of multiplets characterized 
by definite quantum numbers. 

{{Si}B, {Si}M} containing the spacings of, respectively, the baryon, the meson, and all 
multiplets. As usual, the mean spacing, (S) = (1/-/V) Z^^i Sk is scaled out and we deal with 
dimensionless arrays: 

{Si}B = {Si}B/ {S)b , 

{Si}M = {Si}M/ {S)m , (1) 
{SiU = {Si}H/{S)H^{{S,}B/{S)„, {Si}M/{S)H}. 

The mean spacing values {S)^ are computed in Table I. 
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Spectrum 


{S),m MeV 


Baryon 


289.8 ± 18.8 


Meson 


284.3 ± 35.8 


Hadron 


287.1 ± 27.3 



TABLE I: Mean mass spacing. The error bars stem from the PDG quoted errors. 

Histograms of the resulting spacing distributions, p{{si}B), p{{si}M), and p{{si}H) are 
presented in Fig. ||. The curves in the figure show the Poisson-hke distribution [p(s) = 
exp(— s)], and the Wigner surmise [p{s) = f s exp(— They represent two extreme 
regimes: uncorrelated spectrum vs. a strongly correlated one. The figure clearly shows that 
the experimental distribution is of the Wigner-surmise type. We have checked that if one 
does not discriminate the QNs ( "mixed QNs" ) the spacing distribution is Poissonian, so only 
states with the same QNs are correlated. 

The histogram plots represent the distributions in a rather qualitative way because due 
to the low statistics the picture is very sensitive to the choice of the grid. Accounting for 
the error bars makes it even less accurate. It is more instructive to consider integrated 
characteristics, such as the moments of distribution. For our empirical distribution the n-th 
moment is simply calculated as 

1 ^ 

(^")exp = T7E^r, (2) 

i=l 

which obviously is not sensitive to the analysis artifacts such as the choice of energy grid. 
At the same time let us consider a more general form of the Wigner surmise: 

Pw(s) = A/exp(-Ss2) (3) 

where /3 is a parameter — the Dyson index; A and B are constants fixed by the normalization 
conditions: {s^)^ = 1 = (s^)-y^, where, naturally, {s"')^^ = ds s"^ py/{s). Depending on 
the value of (3 expression @ approximates the spacing distributions of the RMT for various 
types of Gaussian ensembles. The most common values are /? = 1, 2, and 4 which represent 
the orthogonal (GOE), unitary (GUE), and symplectic (GSE) ensembles, respectively. 

Computing the first ten, or so, moments of all the relevant distributions we obtain Fig- 
ure 1^. The moments of the Poisson and Wigner distributions (shown by lines) are propor- 
tional to the Gamma-function [e.g., (s")p = T{n + 1), and = {2/y/n)"'T{^n + 1), for 
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FIG. 2: Histograms of the nearest-neighbor mass spacing distribution for hadron states with same 
QNs. Curves represent the Poisson (dashed) and Wigner (soHd) distributions. 
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(3 = 1]. Therefore they all increase rapidly with n and we have used logarithmic scale to 
plot them. 

Data points in the first three panels correspond to the empirical distributions shown by 
the histograms in Fig. 0, now with the PDG error bars properly accounted for|TP|. The 
Wigner surmise with (3=1 fits these distributions extremely well: comparing the first 
ten moments, x^/Y>oini values for the three distributions are shown in the "None" row of 
Table |n|. 

In general, before combining the level statistics from two different spectra we would need 
to perform an unfolding of these spectra. Then the combined meson-baryon array of spacings 
could be constructed as {si}'^ = {{si}B, {sjIm}- liowever, because the baryon and meson 
mean spacing come out so similar (Table I) this construction leads to practically same results 
as the one in Eq. (|^), the value changes insignificantly, see "B vs M" row of Table ||. 
One can also perform the unfolding with respect to strangeness, which means one treats 
the spectrum of strange baryons and mesons as different from that of non-strange ones, 
and unfolds {i.e., normalizes to unit mean spacing) prior combining them. In this case the 
agreement with Wigner surmise has slightly improved as is seen from the "S=0 vs —1" row 
of Table 0. 



Unfolding 


Baryon 


Meson 


Hadron 


None 


0.52 


0.66 


0.57 


B vs M 






0.54 


S=0 vs -1 


0.29 


0.40 


0.33 



TABLE II: /point fit of /3 = 1 Wigner surmise to the first ten moments of the experimental 
spacing distributions for various unfoldings. 



To recapitulate, the hadron level-spacing distribution is remarkably well described by 
the Wigner surmise for (3 = 1. This indicates that the fluctuation properties of the hadron 
spectrum fall into the GOE universality class, and hence hadrons exhibit the quantum chaos 
phenomenon. One then should be able to describe the statistical properties of hadron 
spectra using RMT with random Hamiltonians from GOE that are characterized by good 
time-reversal and rotational symmetry. 
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The hadron-exchange models (e.5f.[|T^) should also be able benefit from the RMT meth- 
ods, in a way analogous to the stochastic theory of compound-nucleus reactions |18 . 
However, since unitarity is an essential ingredient of such models, a study of hadron width 
distribution is in order. Widths of compound nuclei have the Porter-Thomas distribution 
0] that eventually allows one to define the unitary S'-matrix completely within the RMT. 
We anticipate that the same approach will be applicable in the hadron-exchange reactions 
once the hadron-width distribution is determined. 

It would be interesting to see how quantum chaos emerges in hadrons as a result of quark- 
gluon underlying dynamics. Lattice QCD studies in this direction are underway 0, ||, 
This problem could also be addressed in a quark-model framework. In particular, does a 
quark-model Hamiltonians which fits the physical spectrum indeed defines a chaotic classical 
system? 

In conclusion, we have examined the nearest-neighbor level-spacing (mass-splitting) dis- 
tribution, p{s), of the experimental hadron spectrum. We focused on the lighter part of the 
hadron spectrum (m < 2.5GeV), since it is more reliably known. The mean level-spacing 
seems to be the only relevant scale here, and after it is scaled out the distribution shows 
universal behavior. Unfortunately, the masses and quantum numbers are well known only 
for a few dozen hadron states, so achieving high statistics is out of the question. Nonethe- 
less, the low-statistics analysis we have performed pinpoints the moments of the distribution 
accurately enough to claim that hadronic p{s) fits the Wigner surmise with linear level re- 
pulsion (/3 = 1). This indicates that the spectrum falls into the GOE universality class of 
random matrix theory. Invoking the BGS conjecture, this result is viewed as an empirical 
evidence of the "quantum chaos" phenomenon in hadrons. 
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